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Abstract 

Integrated  quantum  photonics  relies  critically  on  the  purity,  scalability,  integrability,  and  flexibility  of  a 
photon  source  to  support  diverse  quantum  functionalities  on  a  single  chip.  Up  to  date,  it  remains  an  open 
challenge  to  realize  an  efficient  monolithic  photon-pair  source  for  on-chip  application.  Here  we  report 
a  device  on  the  silicon-on-insulator  platform  that  utilizes  dramatic  cavity  enhanced  four-wave  mixing  in 
a  high-Q  silicon  microdisk  resonator.  The  device  is  able  to  produce  high-purity  photon  pairs  in  a  comb 
fashion,  with  an  unprecedented  spectral  brightness  of  6.24  x  107  pair/ s/mW2 /GHz  and  photon-pair  cor¬ 
relation  with  a  coincidence-to-accidental  ratio  of  1386  ±278  while  pumped  with  a  continuous-wave  laser. 
The  superior  performance,  together  with  the  structural  compactness  and  CMOS  compatibility,  opens  up  a 
great  avenue  towards  quantum  silicon  photonics  with  unprecedented  capability  of  multi-channel  parallel 
information  processing  for  both  integrated  quantum  computing  and  long-haul  quantum  communication. 
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Photon-based  quantum  information  technology  has  found  broad  applications  ranging  from 
quantum  communication1,  quantum  computing2,  to  quantum  metrology3.  Recent  advances  in 
integrated  quantum  photonics4"7  show  great  promise  for  chip-scale  quantum  information  process¬ 
ing  with  tremendous  complexity.  A  bright,  single-mode,  high-purity,  integrated  source  of  single 
photons  and/or  entangled  photon  pairs  is  essential  for  all  these  applications8-10,  particularly  for  in¬ 
tegrated  quantum  photonic  circuits  which  rely  critically  on  the  purity,  scalability,  integrability,  and 
flexibility  of  the  photon  source  to  support  diverse  quantum  functionalities  on  a  single  chip11-13. 
However,  although  a  variety  of  photon  sources  have  been  developed  in  the  past  few  decades8-13,  it 
remains  an  open  challenge  to  realize  an  efficient  monolithic  source  for  on-chip  application,  lack  of 
which  current  quantum  photonic  devices  have  to  rely  on  external  sources  for  proper  operation4-7. 
Here  we  propose  and  demonstrate  an  ultra-bright  high-purity  chip-scale  photon-pair  source  on 
the  silicon-on-insulator  (SOI)  platform.  By  taking  advantage  of  the  dramatic  cavity  enhanced 
four-wave  mixing  in  a  high-quality  silicon  microdisk  resonator,  we  are  able  to  achieve  a  spec¬ 
tral  brightness  of  6.24  x  107  pair /s/mW2 /GHz  orders  of  magnitude  larger  than  other  photon-pair 
sources11,14-35,  and  an  unprecedented  photon-pair  correlation  with  a  coincidence-to-accidental 
ratio  of  1386  ±278  while  pumping  with  a  continuous- wave  laser.  In  particular,  the  unique  de¬ 
vice  characteristics  enable  photon  pair  generation  in  multiple  frequency  combs,  thus  significantly 
extending  the  wavelength  management  capability  for  integrated  quantum  photonics.  The  unprece¬ 
dented  device  performance  together  with  its  CMOS  compatibility  now  opens  the  door  towards 
quantum  silicon  photonics  with  great  potential  for  multi-channel  parallel  operation  of  novel  quan¬ 
tum  functionalities  on  chip. 

To  date,  nearly  all  photon-pair  sources  are  based  upon  spontaneous  parametric  down  conver¬ 
sion  (SPDC)  or  four- wave  mixing  (FWM)  in  nonlinear  optical  crystals/waveguides8-11.  Bulk 
crystals  emit  photon  pairs  into  a  multimode  spatial  profile,  resulting  in  a  fairly  low  photon  gen¬ 
eration/collection  efficiency10,11,14,15.  Significant  efforts  have  been  devoted  in  recent  years  to 
developing  waveguide  sources  for  single-mode  emission11,16-21.  However,  the  produced  photon 
pairs  generally  exhibit  a  non-factorable  spectrum  which  degrades  considerably  the  quantum-state 
purity8.  FWM  in  silica  optical  fibers  recently  appears  as  a  promising  approach  with  a  great  flexibil¬ 
ity  of  engineering  photon  spectrum22-27,  which,  unfortunately,  suffers  seriously  from  broadband 
Raman  scattering  of  silica25,26,36. 

Silicon  exhibits  a  strong  Kerr  nonlinearity  for  nonlinear  optical  interaction  and  a  large  refractive 
index  enabling  tight  mode  confinement,  which  has  been  explored  intensively  recently  for  a  vari- 
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FIG.  1 :  a,  Schematic  of  generating  a  photon-pair  comb  from  a  silicon  microdisk  resonator  sitting  on  a  silica  pedestal. 
All  photon  modes  couple  to  the  same  propagation  mode  inside  the  delivery  coupling  waveguide,  which  is  a  single¬ 
mode  tapered  silica  fiber,  b,  Group- velocity  dispersion  of  vertically  polarized  transverse-magnetic-like  (TM-like) 
cavity  modes  with  four  different  radial  orders,  termed  as  Po,  Pi,  P2,  and  P3,  respectively,  for  a  silicon  microdisk  with 
a  thickness  of  260  nm  and  a  radius  of  4.9  /im,  simulated  by  the  finite  element  method  (Appendix  A).  The  insets  show 
the  simulated  optical  field  profiles. 

ety  of  applications37-40.  In  particular,  single-crystalline  silicon  has  a  clean  phonon  spectrum  with 
Brillouin-zone-center  phonons  of  a  well-defined  frequency  of  15.6  THz  and  a  narrow  linewidth  of 
~105  GHz,  which  eliminates  the  deleterious  broadband  Raman  noises37,41.  These  superior  fea¬ 
tures  together  with  mature  nanofabrication  technology  make  SOI  an  ideal  platform  for  integrated 
quantum  photonic  application28-35. 

The  device  structure  we  employ  for  high-quality  photon-pair  generation  is  a  compact  silicon 
microdisk  resonator  sitting  on  a  silica  pedestal  (Fig.  la,  with  the  scanning-electron  microscopic 
image  shown  in  Fig.  2).  An  excellent  feature  of  a  microdisk  is  that  its  group-velocity  disper¬ 
sion  is  dominantly  determined  by  the  disk  thickness,  allowing  flexible  dispersion  engineering  for 
the  FWM  process.  As  shown  in  Fig.  lb,  with  a  disk  thickness  of  260  nm,  the  zero-dispersion 
wavelength  (ZDWL)  for  the  vertically  polarized  transverse-magnetic-like  (TM-like)  mode  can  be 
tailored  to  the  useful  telecommunication  band  around  1.5  /im.  Detailed  analysis  shows  that  the 
dispersion  characteristics  support  phase  matching  over  a  broad  spectrum  around  100  nm,  thus  en¬ 
abling  simultaneous  generation  of  photon  pairs  in  a  comb  fashion  over  all  phase-matched  cavity 
modes. 
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In  particular,  the  microdisk  resonator  uniquely  supports  multiple  mode  families  with  different 
radial  orders  (Fig.  lb,  inset).  As  they  experience  different  mode  confinement  introduced  by  the 
whispering-gallery  geometry,  these  mode  families  exhibit  ZDWLs  shifted  by  about  30-40  nm  from 
each  other  (Fig.  lb).  Consequently,  a  variety  of  desired  photon-pair  combs  can  be  produced  inside 
a  single  device  simply  by  pumping  at  different  mode  families.  Moreover,  the  complete  air  cladding 
ensures  a  clean  interaction  of  the  optical  field  with  the  silicon  core  and  thus  eliminates  any  potential 
contamination  from  the  cladding  material,  in  contrast  to  waveguides/microrings28-3 1,33-35  where 
the  optical  interaction  with  the  buried  oxide  layer  is  likely  to  produce  Raman  noise  photons. 

Most  importantly,  the  high-Q  cavity  drastically  modifies  the  density  of  states  of  vacuum,  re¬ 
sulting  in  significant  Purcell  enhancement  on  cavity  modes  which  not  only  dramatically  increases 
the  generation  efficiency  of  photon  pairs  but  also  ensures  them  to  be  created  at  discrete  single 
frequencies  with  a  high  state  purity42.  With  Np  pump  photons  inside  the  cavity,  the  probability  of 
emitting  a  pair  of  signal  and  idler  photons  at  time  ts  and  //  into  the  coupling  waveguide  (Fig.  la)  is 
given  by  (Appendix  B) 

PcUs-ti)  =  (gNp)2e~r' j (1) 


with  rtj  =  rts  (or  rti)  on  the  exponent  when  ts  >  ti  (or  ts  <  hi  where  rej  and  rtJ  (j  =  s.i)  are  the 
external  coupling  rate  and  photon  decay  rate  of  the  loaded  cavity,  for  the  signal  and  idler  photons, 
respectively,  and  T  =  (rts  +  rtj)/2  represents  the  average,  g  =  describes  the  vacuum 

coupling  rate  of  the  FWM  process,  where  ri2,  V,  and  t]  are  the  Kerr  nonlinear  coefficient,  the 
effective  mode  volume,  and  the  spatial  mode  overlap  among  the  interacting  photons,  respectively. 

As  indicated  by  Eq.  (1),  the  signal  and  idler  photons  are  created  in  a  time  scale  of  the  cavity 
photon  lifetime  within  which  the  photon  pair  remains  highly  correlated.  The  longer  the  photon 
lifetime,  the  larger  the  pair  emission  probability.  The  resulting  photon-pair  emission  rate  is  given 
by 


Rc 


resrei2  (gNpf 
rtsrti  r 


(2) 


FWM  creates  photon  pairs  inside  the  cavity  at  a  rate  of  2(gNp)2 /T,  which  are  delivered  to  the 
coupling  waveguide  with  a  pair  extraction  efficiency  of  resrei/ (rtsrti).  A  high  cavity  Q,  a  strong 
Kerr  nonlinearity,  and  a  small  effective  mode  volume  would  result  in  a  large  photon-pair  emission 
rate.  For  example,  simulations  by  the  finite-element  method  show  that  a  silicon  microdisk  with 
a  radius  of  5  gm  exhibits  a  small  effective  mode  volume  of  ~  7  gm3  for  a  Pq  mode,  resulting 
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FIG.  2:  Schematic  of  the  experimental  setup  for  optical  spectroscopy  and  coincidence  photon  counting.  The  pump 
laser  is  coupled  into  the  device  via  a  single-mode  tapered  silica  fiber.  Before  coupling  into  the  device,  the  pump  laser 
passes  through  a  bandpass  filter  and  a  coarse  wavelength-division  multiplexer  (CWDM  MUX)  to  cut  the  laser  noises. 
The  CWDM  MUX  has  a  3-dB  bandwidth  of  17  nm  for  each  of  its  transmission  bands  whose  center  wavelengths  are 
separated  by  20  nm  apart  with  a  band  isolation  of  >  120  dB.  The  bandpass  filter  is  identical  to  one  band  of  the  CWDM 
MUX.  The  produced  photon-pair  comb  is  separated  into  individual  photon  modes  by  the  CWDM  DEMUX  which  is 
identical  to  the  CWDM  MUX  used  at  the  input  end.  The  photoluminescence  spectrum  of  the  photon-pair  comb  is 
recorded  at  each  transmission  port  of  the  CWDM  DEMUX  for  easy  suppression  of  the  pump  wave.  For  coincidence 
counting,  the  photon  pairs  are  recorded  by  gated  InGaAs  single  photon  detectors  (SPDs).  Two  tunable  bandpass 
filters  with  a  3-dB  bandwidth  of  1.2  nm  are  used  to  cut  the  Raman  noises  produced  by  the  delivery  silica  fiber.  The 
pump  power  is  controlled  by  variable  optical  attenuators  (VOAs)  and  its  wavelength  is  calibrated  by  a  Mach-Zehnder 
interferometer  (MZI).  The  inset  shows  the  scanning  electron  microscopic  image  of  the  device. 

in  a  FWM  coupling  rate  of  g  ~  2n  x  460  Hz.  Therefore,  even  a  small  pump  power  of  10  /iW 
launched  into  a  critically  coupled  cavity  with  a  Q  factor  of  106  is  able  to  emit  photon  pairs  at  a 
rate  ^6x  106  pair/s,  clearly  showing  the  extremely  high  efficiency  of  the  proposed  device. 

The  fabricated  device  (see  Method)  was  tested  using  the  experimental  setup  shown  in  Fig.  2. 
For  easy  separation  of  individual  photon  modes  within  the  generated  photon  comb,  the  device  was 
designed  with  a  free  spectral  range  (FSR)  close  to  20  nm  such  that  the  cavity  resonances  fit  directly 
into  the  transmission  bands  of  a  standard  coarse  wavelength-division  multiplexer  (CWDM).  The 
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FIG.  3:  a,  The  transmission  spectrum  of  the  passive  cavity  scanned  by  two  tunable  lasers  operating  at  different 
spectral  ranges  (indicated  as  blue  and  green).  The  two  sets  of  arrows  indicate  cavity  modes  with  constant  FSRs  within 
the  Pi  and  P2  mode  families,  respectively.  The  insets  show  detailed  transmission  spectra  of  a  P2  mode  at  1497.2  nm 
and  a  Pi  mode  at  1564.9  nm,  with  theoretical  fitting  shown  in  red.  b,c,  The  photoluminescence  spectra  of  the  P2 
and  Pi  photon-pair  comb,  respectively,  with  the  pump  nearly  critically  coupled  to  the  cavity.  The  colors  indicate  the 
individual  spectra  recorded  at  different  transmission  ports  of  the  CWDM  DEMUX.  The  four  mode  pairs  in  the  two 
combs  are  termed  as  P2,  P2,  P’ ,  and  P1/,  respectively,  as  indicated  on  the  figures.  The  spectrum  of  the  Pi  comb  extends 
to  the  photon  mode  at  1584.1  nm  which  is  beyond  our  laser  scanning  range  and  is  thus  not  seen  in  a. 

pump  wave  and  the  created  photon  pairs  are  delivered  into  and  out  of  the  microdisk  through  a 
tapered  fiber  which  provides  flexible  control  of  light  coupling. 

The  laser-scanned  transmission  spectrum  of  the  passive  cavity  (Fig.  3a)  shows  four  clear  high¬ 
er  mode  families.  A  constant  FSR  of  2.4225  THz  (corresponding  to  ~18  nm  around  1.5  gm) 
is  observed  to  maintain  among  five  cavity  modes  for  the  P2  mode  family,  covering  the  spectrum 
from  1478  to  1555  nm.  Similar  phenomenon  is  observed  for  the  Pi  mode  family  with  a  constant 
FSR  of  2.3158  THz,  over  the  spectral  region  beyond  1510  nm.  A  detailed  scan  of  the  P2  mode  at 
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1497.2  nm  and  the  Pi  mode  at  1564.9  nm  (Fig.  3a,  insets)  shows  intrinsic  Q  factors  of  3.47  x  105 
and  4.94  x  105,  respectively.  Cavity  modes  within  a  same  mode  family  exhibit  similar  optical 
Q  since,  with  a  same  mode  profile,  they  are  of  similar  sensitivity  to  scattering  loss  induced  by 
fabrication  imperfection.  These  two  mode  families  are  used  below  for  photon-pair  generation. 
The  Po  and  P3  mode  families  are  not  used  in  experiment  since  their  ZDWLs  (~1572  and  1466  nm, 
Fig.  lb)  are  outside  the  scanning  ranges  of  our  lasers,  although  the  Po  mode  family  exhibits  a 
higher  intrinsic  Q  factor  of  about  1  x  106. 

The  constant  mode  spacing  together  with  a  high  optical  Q  ensures  efficient  generation  of  a 
comb  of  photon  pairs.  As  shown  by  the  photoluminescence  (PL)  spectra  in  Fig.  3b,  pumping  at 
the  P2  mode  located  at  1515.6  nm  is  able  to  produce  a  comb  of  four  clean  photon  modes  within 
the  P2  mode  family  (Fig.  3b).  Moreover,  switching  to  pump  at  the  Pi  mode  at  1546.3  nm  is 
able  to  produce  another  set  of  clean  photon  comb  within  the  Pi  mode  family  (Fig.  3c).  For  the 
simplicity  of  notation,  we  term  the  photon  mode  pairs  in  the  two  combs  as  Pj,  Pj1,  Pi,,  and  P^1, 
respectively  (Fig.  3b, c).  The  spectrum  of  each  photon  mode  is  so  sharp  that  it  is  beyond  the 
resolution  of  our  spectrometer  (~  0.135  nm),  implying  the  high  coherence  of  generated  photons. 
The  amplitude  difference  within  a  photon  comb  is  primarily  due  to  different  external  couplings 
of  cavity  modes  to  the  tapered  fiber.  When  the  pump  mode  is  critically  coupled  to  the  cavity, 
the  signals  at  shorter  wavelengths  are  under  coupled  while  the  idlers  at  longer  wavelengths  are 
over  coupled,  resulting  in  a  higher  photon  extraction  efficiency  for  the  idler  within  each  mode  pair 
(Fig.  3b, c).  It  also  leads  to  a  smaller  overall  loaded  Q  for  the  PJ,1  (P1/)  mode  pair  compared  with 
P2  (Pi),  thus  resulting  in  a  lower  pair  production  rate.  Note  that  the  PL  spectra  also  show  some 
other  modes  with  small  amplitudes  whose  physical  origin  is  not  clear  at  this  moment.  Potential 
origins  are  multiphonon  scattering  in  silicon43  or  light  scattering  from  certain  defects  on  the  device 
surface,  which  require  further  investigation.  As  these  tiny  modes  are  well  separated  from  the 
FWM-created  photon  combs,  they  have  negligible  effect  on  the  comb  performance. 

To  characterize  the  performance  of  the  photon-pair  combs,  we  carried  out  coincidence  photon 
counting  for  each  mode  pair  when  the  device  is  pumped  with  a  continuous-wave  (CW)  laser.  As 
shown  by  Fig.  4e-h,  the  coincidence  counts  for  all  correlated  mode  pairs  depend  quadratically  on 
the  pump  power,  indicating  they  are  all  produced  through  the  FWM  process.  The  experimentally 
recorded  pair  fluxes  (blue  dots)  are  close  to  the  theoretical  predictions  (red  curves)  (Appendix  B), 
with  a  small  discrepancy  primarily  coming  from  the  uncertainty  of  Kerr  nonlinear  coefficient  as 
well  as  measured  optical  Q.  No  significant  saturation  is  observed,  thus  indicating  the  negligible  im- 
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FIG.  4:  Recorded  photon  emission  fluxes,  pair  emission  fluxes,  and  CARs  for  the  photon-pair  combs.  Each  column 
of  figures  corresponds  to  one  mode  pair  as  indicated  on  the  top.  The  figures  within  each  row  are  plotted  with  identical 
vertical  axes  for  each  comparison  among  different  mode  pairs,  a-d,  Recorded  emission  fluxes  of  individual  photon 
modes  within  each  mode  pair,  with  detector  dark  counts  subtracted.  The  signal  (idler)  at  the  shorter  (longer)  wave¬ 
length  is  shown  as  blue  (red).  Solid  curves  show  the  theoretical  prediction  of  pure  FWM-created  photon  fluxes,  e-h, 
Recorded  emission  fluxes  of  correlated  photon  pairs.  In  each  figure,  the  left  and  right  vertical  axes  show  the  emission 
flux  and  the  corresponding  true  coincidence  count  per  gate,  respectively,  indicated  by  brown  and  green  colors.  The 
solid  curves  show  the  theoretical  prediction,  i-1,  Recorded  CAR  for  each  mode  pair.  No  any  accidental  coincidence 
is  subtracted.  The  red  solid  curves  show  the  theoretical  prediction.  The  black  dashed  curves  are  calculated  from 
independently  recorded  true  coincidence  counts  and  photon  fluxes  of  individual  modes  (Appendix  C).  All  data  were 
recorded  with  a  detector  gating  width  of  10  ns,  a  quantum  efficiency  of  15%,  a  clock  frequency  of  250  kHz,  and  a 
dead  time  equal  to  the  clock  period. 

pact  of  free-carrier  absorption  and  two-photon  absorption.  In  particular,  a  small  input  pump  power 
of  79  /I W  is  able  to  generate  a  photon-pair  flux  of  8.55  x  105  pair/s  at  PJ>  (Fig.  4e),  clearly  show¬ 
ing  the  high  efficiency  of  the  device.  A  metric  characterizing  the  photon-pair  generation  efficiency 
is  the  spectral  brightness ,  defined  as  the  pair  flux  per  unit  spectral  width  per  unit  pump  power 
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square  (since  FWM  depends  quadratically  on  the  pump  power).  The  spectral  width  of  each  photon 
mode  is  determined  by  the  linewidth  of  the  loaded  cavity,  which  is  1.59  and  2.80  GHz  for  the  PJ> 
mode  pair  at  1497.2  and  1534.3  nm,  respectively.  It  infers  that  the  spectral  brightness  of  the  gen¬ 
erated  photon  pairs  is  6.24  x  107  pair/s/mW2/GHz,  which  is  about  3  orders  of  magnitude  larger 
than  any  FWM-based  photon-pair  sources  reported  up  to  date21-35.  P" ,  P1 ,  and  P'1  exhibit  similar 
spectral  brightness,  with  magnitudes  of  1.78  x  107,  2.32  x  107,  and  2.57  x  107  pair/s/mW2/GHz, 
respectively. 

Indeed,  the  cavity  enhancement  is  so  strong  inside  the  device  that  the  demonstrated  photon- 
pair  source  is  even  brighter  than  SPDC-based  sources,  although  FWM,  a  third-order  nonlinear 
process,  is  generally  much  weaker  than  SPDC  which  is  a  second-order  nonlinear  process.  The 
spectral  brightness  of  SPDC-based  sources  is  generally  described  by  the  slope  efficiency  given  as 
the  photon-pair  flux  per  unit  spectral  width  per  unit  power15-18.  P|>  shown  in  Fig.  4e  has  a  slope 
efficiency  of  4.93  x  106  pair/s/mW/GHz  at  a  pair  flux  of  8.55  x  105  pair/s,  which  is  more  than 
one  order  of  magnitude  larger  than  any  SPDC-based  waveguide  sources11,14-20.  It  is  comparable 
with  the  state-of-the-art  cavity-enhanced  SPDC  sources44^16,  while  the  latter  create  photon  pairs  in 
multi-spatiotemporal  modes  which  require  stringent  filtering.  Moreover,  the  latter  are  constructed 
with  bulk  optical  components  which  require  sophisticated  cavity  locking  and  are  challenging  for 
chip-scale  integration. 

Figure  4  shows  that  the  pair  fluxes  are  smaller  than  those  of  corresponding  individual  photon 
modes.  This  is  because  of  the  non-unity  photon  extraction  efficiency,  a  common  feature  of  cavity 
quantum  electrodynamic  systems:  Although  the  signal  and  idler  photons  are  always  created  in  a 
pair  fashion  inside  the  cavity,  they  transmit  out  of  the  cavity  independently,  only  a  certain  fraction 
of  which  remains  the  pair  correlation.  The  performance  can  be  improved  by  using  cavity  modes 
with  higher  intrinsic  Q  ( e.g .,  the  Po  mode  family)  together  with  higher  external  coupling.  On  the 
other  hand,  Fig.  4a-d  show  that  the  fluxes  of  individual  photon  modes  deviate  slightly  away  from 
the  quadratic  power  dependence,  indicating  that  they  are  accompanied  with  a  certain  amount  of 
noises.  Detailed  characterizations  show  that  the  noise  photons  dominantly  come  from  the  broad¬ 
band  Raman  noises  generated  inside  the  delivery  fibers  (Appendix  D).  In  experiment,  a  pair  of 
bandpass  filters  with  a  3-dB  bandwidth  of  1.2  nm  were  used  in  front  of  the  single  photon  detectors 
(SPDs)  (Fig.  2)  to  suppress  the  Raman  noises.  As  each  photon  mode  has  a  very  narrow  spectral 
width,  the  Raman  noises  can  be  suppressed  further  by  using  a  narrower  filter.  They  can  also  be 
reduced  by  shortening  the  lengths  of  delivery  fibers.  Such  Raman  noises  would  be  absent  if  an 
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on-chip  laser39  is  used  to  generate  the  photon-pair  combs. 

An  important  figure  of  merit  characterizing  the  quantum  correlation  between  the  photon  pairs 
is  the  coincidence-to-accidental  ratio  (CAR)18,24,26,31,  which  is  directly  related  to  the  quantum- 
interference  visibility  of  the  entangled  states  constructed  from  the  photon  pairs36.  Figure  4i-l 
show  the  recorded  CAR  (blue  dots)  for  all  mode  pairs.  In  general,  CAR  increases  with  decreased 
photon-pair  flux  because  of  reduced  probability  of  multi-pair  creation  and  saturates  at  a  certain 
level  of  pair  flux  when  the  detector  dark  counts  start  to  dominate.  All  mode  pairs  exhibit  very  high 
CAR  at  all  power  levels.  For  the  P2  mode  pair,  the  device  is  able  to  achieve  unprecedented  CAR 
with  a  value  of  274  ±  10  at  a  pair  flux  of  8.55  x  105  pair/ s  and  a  peak  value  of  1386  ±278  at  a 
pair  flux  of  2.35  x  104  pair/s.  These  values  are  more  than  three  times  larger  than  other  photon- 
pair  sources17-35,  even  without  using  pulsed  pump  to  suppress  noise  photons18,24,25,28,29,31,32  or 
delicate  superconducting  SPDs  to  reduce  detector  dark  counts16,26,34,35.  Similar  CAR  is  achieved 
for  P1,  with  values  of  1097  ±  165  and  244±6  at  pair  fluxes  of  2.43  x  104  and  6.31  x  105  pair/s, 
respectively.  At  the  same  time,  P11  and  P^1  exhibit  as  well  relatively  high  CAR,  with  peak  values  of 
733  ±  106  and  679  ±74,  respectively.  For  all  the  mode  pairs,  the  power-dependent  trend  of  CAR 
follows  closely  the  theory  (red  curves)  (Appendix  C).  The  discrepancy  primarily  comes  from  the 
uncertainty  in  characterizing  Raman  noises  generated  in  delivery  fibers.  These  results  clearly  show 
the  superior  quality  of  generated  photon-pair  combs.  We  expect  that  the  peak  CAR  of  the  photon- 
pair  combs  can  be  improved  by  at  least  one  order  of  magnitude  in  the  future,  by  suppressing  the 
Raman  noises  produced  in  the  delivery  fibers  and  by  using  better  detectors  ( e.g .,  superconducting 
SPDs)  with  lower  dark  counts. 

The  high  pair  correlation  and  high  coherence  of  the  photon-pair  combs  allow  resolving  the 
temporal  structure  of  quantum  correlation  between  the  signal  and  idler  photons.  This  is  particularly 
enabled  by  the  tapered-fiber  coupling  system  which  provides  flexible  control  of  the  external  light 
coupling  and  thus  the  cavity  photon  lifetime.  Figure  5  shows  the  normalized  coincidence  spectra 
of  P3>  with  two  different  external  coupling  conditions.  When  the  tapered  fiber  is  close  to  the  device, 
the  external  coupling  is  strong  which  results  in  a  short  photon  lifetime  of  65  ps  (averaged  between 
signal  and  idler).  As  a  result,  the  coincidence  spectrum  (blue  curve)  shows  a  full  width  at  half 
maximum  (FWHM)  of  370  ps,  which  is  19%  larger  than  the  time  resolution  of  the  instrument 
response  function  (IRF,  gray  curve)  of  our  coincidence  counting  system.  However,  moving  the 
tapered  fiber  away  from  the  device  increases  the  photon  lifetime  considerably  to  171  ps,  leading 
to  a  coincidence  spectrum  (green  curve)  with  a  FWHM  of  453  ps  which  is  45%  larger  than  the 
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FIG.  5:  Normalized  delay-time-dependent  coincidence  spectra  for  Pl2.  The  blue  and  green  curves  show  the  cases 
with  different  external  coupling  conditions  and  thus  different  averaged  photon  lifetimes.  The  gray  curve  shows  the 
instrument  response  function  (IRF)  of  our  coincidence  counting  system,  with  a  FWHM  of  312  ps  which  primarily 
come  from  the  timing  jitters  of  the  two  SPDs.  The  red  and  purple  curves  show  the  theoretical  prediction.  The 
coincidence  spectra  were  recorded  with  a  detector  gating  width  of  10  ns,  a  quantum  efficiency  of  25%,  and  a  clock 
frequency  of  125  kHz.  The  coincidence  counter  has  a  time  resolution  of  4  ps. 

time  resolution  of  IRF.  The  experimentally  observed  time  correlation  is  explained  well  by  the 
theory  (Appendix  C).  To  the  best  of  our  knowledge,  this  is  the  first  time  to  resolve  the  temporal 
photon  correlation  for  a  chip-scale  photon-pair  source,  although  it  was  observed  before  in  cavity- 
enhanced  SPDC  sources44-46.  Quantum  coherence  between  correlated  photon  pairs  is  essential  for 
constructing  a  variety  of  quantum  functionalities1,2’9.  The  photon  lifetime  of  171  ps  in  our  device 
corresponds  to  a  propagation  length  of  ~  1 .5  cm,  which  is  about  the  footprint  of  a  photonic  chip. 
Therefore,  the  produced  photon  pairs  are  able  to  share  temporal  coherence  over  the  entire  chip, 
which  would  significantly  improve  the  quality  of  quantum  functionalities  built  upon. 

The  superior  performance  of  the  demonstrated  device  offers  a  great  opportunity  for  realizing 
a  variety  of  high-quality  quantum  functionalities  on  chip.  Straightforward  applications  include 
producing  heralded  single  photons10  and  quantum  entanglement  in  various  degrees  of  freedom 
(frequency47,  path9,  time1,  polarization31,  etc.),  which  are  essential  for  quantum  communication 
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and  quantum  computing1,2.  Compared  with  conventional  approaches,  the  high  photon  coherence 
in  the  device  would  significantly  improve  the  purity  of  the  constructed  quantum  states.  As  the 
cavity  enhancement  only  occurs  within  the  cavity  photon  lifetime,  the  demonstrated  approach  is 
able  to  operate  with  a  repetition  rate  above  gigahertz  by  use  of  pulsed  pumping.  Further  increase 
of  photon-pair  coherence  and  spectral  brightness  can  be  obtained  by  improving  optical  Q  ( e.g ., 
through  surface  passivation48).  Such  a  high  coherence  would  enable  teleporting  quantum  entan¬ 
glement  intra-chip,  inter-chip,  or  off-chip  over  far  distance  free  from  the  stringent  synchronization 
requirement  in  conventional  schemes1,9,49  (The  direct  photon  coupling  into  optical  fibers  with 
negligible  losses  (Fig.  2)  is  particularly  suitable  for  this  application).  Moreover,  the  comb  nature 
of  emitted  photons  makes  readily  available  wavelength-division  multiplexing  in  quantum  regime, 
and  thus  significantly  extends  the  spectral  efficiency  and  wavelength  management  capability  of 
integrated  quantum  photonics.  In  particular,  it  exhibits  great  potential  for  realizing  multipartite 
entanglement9  and  cluster  state50  inside  a  single  device,  which  would  dramatically  enhance  the 
capability  of  integrated  quantum  computing.  On  the  other  hand,  the  high  coherence  and  spectral 
flexibility  of  the  demonstrated  approach  enable  excellent  interfacing  with  other  chip-scale  infor¬ 
mation  processing/storage  elements12.  These  features  together  with  the  developed  silicon  photonic 
technology39,40  offer  a  great  opportunity  to  ultimately  form  a  complete  CMOS  compatible  chip- 
scale  platform  for  quantum  information  processing. 

Methods 

The  silicon  microdisk  resonator  is  fabricated  from  a  standard  silicon-on-insulator  wafer  by  use  of  e-beam  lithogra¬ 
phy  to  define  the  device  pattern,  fluorine-based  plasma  to  etch  the  silicon  layer,  and  hydrofluoric  acid  to  undercut  the 
silica  pedestal.  Light  is  coupled  into  and  out  of  the  microdisk  through  a  tapered  optical  fiber  which  is  anchored  onto 
two  nanoforks  fabricated  near  the  microdisk  for  stable  operation  (see  Fig.  2).  Touching  to  the  nanoforks  introduces 
a  small  insertion  loss  of  only  about  7%,  as  indicated  by  the  cavity  transmission  spectrum  given  in  Fig.  3a.  The  PL 
spectra  were  recorded  at  the  transmission  ports  of  the  CWDM  DEMUX  to  suppress  the  residual  pump  wave  which 
otherwise  would  cause  significant  stray  light  inside  the  spectrometer. 

The  CAR  is  obtained  by  the  following  procedure:  First,  at  each  pump  power,  a  delay-time-dependent  coincidence 
counting  histogram  (we  term  as  coincidence  spectrum  here  for  simplicity.  See  Fig.  5  for  an  example)  is  obtained  with 
zero  time  delay  between  the  two  SPDs.  The  coincidence  spectrum  is  integrated  over  its  full  width  at  half  maximum 
Tfwhm  to  obtain  the  total  coincidence  C.  The  statistical  error  of  C  is  given  by  Ec  =  \fC  since  the  photon  counting 
follows  a  Poisson  distribution.  To  find  the  accidental  coincidence,  a  series  of  coincidence  spectra  are  obtained  with 
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a  time  delay  of  jTo  between  the  two  SPDs,  where  To  is  the  detector  clock  period  and  j  =  1,2,  ...,200.  Integrate  each 
of  them  within  the  same  Tfwhm  to  obtain  an  accidental  coincidence  Aj.  The  average  and  the  standard  deviation  of  Aj 
give  the  accidental  coincidence  A  and  its  statistical  error  E\.  Note  that  no  any  accidental  coincidence  is  subtracted. 
Finally,  the  CAR  is  given  by  CAR  =  (C  —  A)  /A  and  its  statistical  error  Eqar  is  given  by  =  ]J  {^c^j  +  . 

To  reduce  the  statistical  error,  each  data  point  in  Fig.  4  was  recorded  with  a  time  period  of  5400  —  18900  seconds 
depending  on  the  magnitude  of  photon-pair  flux.  With  such  amount  of  photon  counting  data,  ^  is  very  small  (<  2%) 
and  |r^jr  is  dominantly  determined  by  ^  since  the  accidental  coincidence  is  very  small  in  our  device. 

To  find  the  photon-pair  emission  flux  shown  in  Fig.  4e-h,  a  same  procedure  is  carried  out  but  now  the  integration  is 
performed  over  the  entire  time  window  of  coincidence  spectra  (20  ns,  with  a  detector  gating  width  of  10  ns)  to  obtain 
the  total  coincidence  Ct,  the  accidental  coincidence  At,  and  their  statistical  errors  Eq  and  E\t.  The  true  coincidence 
is  thus  given  by  Cj  =Ct—At.  Its  statistical  error  EcT  is  given  by  EcT  =  \J  {Eq)2  +  (£aJ2-  Dividing  Cj  by  the  total 
number  of  detector  gates  within  the  data  collection  time  we  find  the  true  coincidence  counts  per  gate  (which  is  shown 
on  the  right  axes  of  Fig.  4e-h).  The  photon-pair  emission  flux  is  obtained  by  calibrating  the  true  coincidence  counts 
per  gate  with  the  clock  frequency,  duty  cycle,  and  quantum  efficiency  of  the  detectors  as  well  as  the  propagation  losses 
from  the  device  to  the  detectors  for  the  photon  pairs.  The  statistical  error  of  pair  emission  flux  is  obtained  accordingly 
from  Ect  . 

The  IRF  of  our  coincidence  counting  system  is  obtained  with  identical  ultrashort  pulses  launched  onto  the  two 
SPDs.  The  ultrashort  pulses  come  from  a  mode-locked  fiber  laser  with  a  pulse  width  (FWHM)  of  240  fs.  They  are 
broadened  to  about  5  ps  by  the  delivery  fibers  before  they  hit  on  the  SPDs.  The  pulse  width  is  much  shorter  than  the 
timing  jitters  of  the  SPDs,  thus  providing  a  fairly  accurate  approach  to  characterize  the  IRF. 

In  all  photon  counting  measurements,  the  detector  gating  width  is  set  to  10  ns  and  the  deadtime  is  equal  to  the 
clock  period.  The  two  InGaAs  SPDs  have  dark  counts  of  5.09  x  10“5  and  3.18  x  10  5  per  gate,  respectively,  at  a 
quantum  efficiency  of  15%  and  a  clock  frequency  of  250  kHz.  The  detector  dark  counts  increase  to  8.96  x  10-5  and 
6.93  x  10-5  per  gate  at  a  quantum  efficiency  of  25%  and  a  clock  frequency  of  125  kHz.  In  both  cases,  the  after-pulsing 
probability  is  less  than  8%. 
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Appendix  A:  Group- velocity  dispersion  of  a  silicon  microdisk  resonator 

For  each  mode  family,  the  silicon  midrodisk  resonator  is  simulated  by  the  finite-element  method 
to  find  the  cavity  resonances  ©(m)  as  a  function  of  mode  number  m  over  a  broad  spectrum.  As 
the  propagation  constant  of  a  cavity  mode  is  well  approximated  by  k  =  m/R  where  R  is  the  device 
radius,  we  obtain  the  dispersion  relation  k(a>)  for  the  cavity  modes.  Fitting  it  with  a  high-order 
polynomial,  we  obtain  the  group- velocity  dispersion  shown  in  Fig.  lb  of  the  main  text. 

In  general,  the  group- velocity  dispersion  of  a  microdisk  resonator  is  close  to  a  slab  waveguide 
with  the  same  thickness,  which  can  be  used  as  a  rough  guidance  for  searching  the  group-velocity 
dispersion  of  a  microdisk  resonator. 


Appendix  B:  Theory  of  photon-pair  generation  in  a  microdisk  resonator 

In  this  section,  we  provide  the  theory  describing  photon  pair  generation  in  a  silicon  microdisk 
resonator.  We  consider  only  the  Kerr  nonlinear  optical  interaction  and  neglect  two-photon  absorp¬ 
tion  and  free-carrier  effect  because  of  their  negligible  effect  in  our  device.  We  consider  first  the 
case  shown  in  Fig.  6  where  a  pump  wave  launched  into  the  resonator  produces  a  pair  of  signal  and 
idler  photon  which  are  transmitted  out  into  the  coupling  waveguide.  The  Kerr  nonlinear  interaction 
among  the  three  cavity  modes  at  frequencies  <Uoj  (  j  =  p.  s,  i )  can  be  described  by  a  Hamiltonian 
H  =  Hq+Hj,  where  Ho  describes  the  passive  cavity  modes  coupled  to  the  bus  waveguide  and  Hj 
describes  the  Kerr  nonlinear  interaction51  inside  the  cavity,  with  the  following  forms 


Hq=  Yj  IhaojaXaj-h^/T^ 
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Ojbje~i(0jt  +  b^ajeiC0jt 
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where  the  intracavity  field  operator  aj  ( j  =  p.sj)  is  normalized  such  that  a  -cij  represents  the 
photon  number  operator  and  bj  is  the  field  operator  of  the  incoming  wave  at  carrier  frequency 
c Oj  inside  the  coupling  waveguide  normalized  such  that  b^bj  represents  the  operator  of  the  input 
photon  flux.  ( b\bs )  =  (b]bj)  =  0  since  only  the  pump  is  launched  into  the  cavity,  bj  satisfies  the 
commutation  relation  of  [bj(t),b^(t')]  =  8(t  —  t').  rej  is  the  external  coupling  rate  of  the  cavity 
mode  at  rnoj  • 

In  Eq.  (B2),  the  first  and  second  terms  describe  the  self-phase  and  cross-phase  modulation 
(SPM  and  XPM)  from  the  pump  mode,  respectively,  and  the  third  and  fourth  terms  govern  the 
four-wave  mixing  (FWM)  process.  We  have  assumed  the  signal  and  idler  are  much  weaker  than 
the  pump  so  that  we  can  neglect  the  SPM  and  XPM  introduced  by  them.  gp,  gpj  (j  =  s,  i ),  and  gpsj 
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FIG.  6:  Schematic  of  photon  pair  generation  in  a  silicon  microdisk  resonator. 


are  the  vacuum  coupling  rate  for  the  SPM,  XPM,  and  FWM,  respectively.  They  are  given  by  the 
following  general  expression 


£ijkl  — 


3tl  ^(Ooi  ft»0j  ft>0k  ©01  fiijkl 

4£0rt,-n7'ft£rt/\/jjkl 
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where  nv  (v  =  i,j,k,l)  is  the  refractive  index  of  silicon  material  at  COfJv  and  %(3>  is  the  third-order 
nonlinear  susceptibility.  V^i  =  (V/V/V^V/)1^4  is  the  average  effective  mode  volume  and  Vv  (v  = 
i,  j,  k,  l )  is  that  at  individual  frequency  £0ov  given  as52 


Vv  = 


dv£r{v,  COqv)  \Ev(v,  ©)v)  |  j' 
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where  £r(r.  £0ov)  is  the  dielectric  constant  of  material  and  the  subscript  si  in  the  integral  denotes 
the  integral  over  the  silicon  core,  r/^  describes  the  spatial  mode  overlap  given  by52 


„  _  Ldrienerje^y^EIEjElE, 

^jkl  ~  r  _  i  1/4  ’ 

{n^umLidreim4) 


(B5) 


where  £rv  =  £r(r,co qv)  and  Ev  =  Ev ( r ,  tOov )  is  the  electrical  field  profile  for  the  cavity  resonant 
mode  at  toov- 

The  vacuum  coupling  rates  for  various  nonlinear  effects  are  thus  obtained  from  Eq.  (B3)  as 
gp  =  gpppp,  gpj  =  gpjjp  0  =  0.  and  gpsi  =  gspip.  In  our  device,  the  three  cavity  modes  are  of  close 

frequencies  and  similar  mode  profiles,  resulting  in 


gp  ~  gp] 


<?psi  —  g  — 


cr\n2h(0p^(0^0i 

nsmV 


(B6) 
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where  77  and  V  denote  T7spip  and  yspip,  respectively,  and  we  have  used  the  definition  of  Kerr  non¬ 


linear  coefficient37:  ri2  =  f  l  2  at  the  pump  frequency.  Using  Eqs.  (Bl)  and  (B2)  and  counting  in 

4  CsQCTlp 


_  3%(3) 
tip 

the  intrinsic  cavity  loss53,  we  obtain  the  following  equations  of  motion  in  the  Heisenberg  picture 
governing  the  wave  dynamics  inside  the  cavity: 


^  — (  / tx)()p  T'tp/2)a.p  +  igcJpdp  +  /  \/rep/?pc  p  +  i  \JY  p ,  (B7) 

^  =  (-/ft) os  -  rts/2)as  +  2 iga\apas  +  iga]a2  +  i^/f^,bse~l(0st  +  i^/f^sus,  (B8) 
^  =  (-ZftJoi  -  rti/2 )cii  +  2igapapcij  +  iga]a2  +  iy/f^bie~mt  +  iy/f^Ui,  (B9) 


where  rQj  and  rtj  =  Toj  +  re|  ( j  =  p.sj)  are  the  photon  decay  rates  of  the  intrinsic  and  loaded 
cavity,  respectively,  uj  is  the  noise  operator  associated  with  intrinsic  cavity  loss,  which  satisfies 
the  commutation  relation  of  [uj(t),Uj(t')\  =  8(t  —  t').  In  Eqs.  (B7)-(B9),  we  have  neglected  the 
pump  depletion  and  the  signal/idler-induced  SPM  and  XPM.  The  transmitted  field  is  given  by53 


fj  ~  bj  + 


(BIO) 


In  general,  the  pump  mode  can  be  treated  as  a  classical  field  and  Eqs.  (B7)-(B9)  can  be  easily 
solved  to  find  the  solutions  for  as  and  a;,  from  which,  together  with  Eq.  (BIO),  we  find  the  emitted 
photon  fluxes  of  signal  and  idler  as 

I  I  .  /'*—)— 00 

J _  |S(<»)|2<i»,  (Bll) 

T  ’  .1  '  f‘  — )— 00 

{fi{t)fiif))  =  ~^J_oo  l5(w)|2 d(°i  (B12) 

and  the  pair  correlation  as 


Pc(ts,ti)  =  (fj(ti)fhts)fs(ts)fi(ti))  -  (flits) fs(ts)){fi  Wiki)) 


resrei 


1  /*  — )— 00 

—  /  B(-co)  [TtsA(ft))  -  1]  da 

2n  7-oo 


where  A(a>)  and  B(co)  have  the  following  expressions 

.  /  \  _  Tti/2  /O) 

'  =  (r„/2  -  iffl)(rti/2  -  la)  -  (sjVp)2  ’ 

—igci2p 

5(w)  =  (rts/2  -  /ft))(rti/2  -  /©)  -  (g/\g2  ’ 

where  Np  =  (apcip)  is  the  average  photon  number  of  the  pump  wave  inside  the  cavity. 


(B13) 


(B14) 

(B15) 
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Equations  (B1 1)-(B13)  are  quite  general  since  they  include  the  multiphoton  generation  induced 
by  the  stimulated  FWM.  For  single  photon-pair  generation,  Eqs.  (Bll)  and  (B12)  reduce  to 

(fHt)fs(t))  -  2rei(g^)2,  (B16) 

1  ts1 

(fht)fi(t))  ~  2Fe^p)2,  (B17) 

1  ti1 

where  T  =  (rts  +  rti)/2  represents  the  average  photon  decay  rate  of  the  loaded  cavity.  The  photon- 
pair  correlation  reduces  to 


Pc(ts,ti )  ~ 


gNp)2e  r‘s^  ti]  if  ts  >  ti 
{gNp)2e~T^ti~ts">  if  ts  <  ti 


(B18) 


Equation  (B18)  can  be  regarded  as  the  joint  probability  density  of  emitting  a  photon  pair  into  the 
coupling  waveguide.  It  is  given  as  Eq.  (1)  in  the  main  text.  Integrating  Eq.  (B18)  over  T  =  ts  —  ti, 
we  obtain  the  average  photon-pair  emission  flux  as 


Rc 


resrei2  (gNpf 

rtsrti  r 


(B19) 


which  is  Eq.  (2)  in  the  main  text.  Comparing  Eq.  (B19)  with  Eqs.  (B16)  and  (B17),  we  can 
find  that  the  photon-pair  generation  rate  inside  the  cavity  is  2{gNp)2 /T  and  the  photon  extraction 
efficiency  is  rej/rtj  (j  =  s,i)  for  the  individual  signal  and  idler  mode  while  it  is  resrei/(rtsrti) 
for  the  correlated  photon  pair.  As  rej/Ttj  is  less  than  1,  the  photon-pair  flux  is  smaller  than  the 
individual  photon  fluxes.  This  is  experimentally  observed  in  Fig.  4  in  the  main  text. 

In  practice,  the  situation  is  much  more  complicated  than  the  ideal  situation  described  above, 
since,  due  to  the  light  scattering,  the  optical  wave  inside  the  cavity  will  couple  to  the  degenerate 
mode  propagating  backward  in  the  opposite  direction  (Fig.  7).  For  a  passive  cavity,  this  phe¬ 
nomenon  manifests  as  a  doublet  in  the  transmission  spectrum  (Insets  of  Fig.  3a  in  the  main  text). 
In  this  case,  The  Hamiltonian  Hq  and  H[  become 


H0=  {^j  (fljffljf  +  «jbajb)  “  +  bl PJ flj 


,baif 


J=P,S,l 


Ht  = 


~hV^}  (alfbif  +  albbib)e  m,t  +  (^jf^jf  +  blba]b)elb¥ 
k8p  half? ali +  (4b)2flpb  +  4«Jf«pf4bflpb 


2h(a^fapf  ^pb^pb)  §ps(^Jf^sf  T  ^sb^sb )  T  gpi  T  ^\ba\b ) 


-frg 


psi 


t  t  2  i  t  t  2 
asfalfapf  +  asbalbapb 


-  %psi  (aif)2asfaif  +  (alb)2asbaib 


V 


pb' 


(B20) 


(B21) 
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FIG.  7 :  Schematic  of  photon  pair  generation  in  a  silicon  microdisk  resonator.  The  clockwise  and  counter-clockwise 
modes  are  coupled  inside  the  cavity. 


where  ajf  and  cv^  (j  =  p,s,i)  represent  the  optical  fields  propagating  clockwise  and  counter¬ 
clockwise,  respectively,  inside  the  cavity  (Fig.  7),  and  /3 j  is  the  coupling  coefficient  between  them. 
bjf  and  represent  the  input  fields  propagating  forward  (input  from  the  left  end)  and  backward 
(input  from  the  right  end),  respectively,  inside  the  coupling  waveguide.  The  transmitted  fields  now 
become  fj  =  bjf  +  i^/T^ajf.  Using  Eqs.  (B20)  and  (B21)  and  following  the  same  procedure,  we 
can  find  the  equations  of  motion  as 

,  =  (  io Op  —  rtp/2)flpf  +  ifipUpb  +  ig  (rZppCZpf  +  2r/pbflpb  )r/pf  -F  f  Cpf  ?  (B22) 

=  (  Dp  rtp/2)apb  +  ifipCi pf  +  ig(a pbOpb  d-  20pfQpf)opb  d-  ?Cpb •>  (B23) 

=  (-iCQos  ~  rts/2)asf  +  //3v«sh  +  2/g(aJfapf  +  ajbapb)asf  +  iga^a^  +  i’Csf,  (B24) 

^  =  (  i C0qs  rts/2)asb  T  ifis  asf  -F  2ig(ci^cipf  T  Opb£?pb)risb  T  (g£(b^pb  d-  i^sb-  (B25) 

—  =  (-10)01  -rti/2)flif  +  iA«ib  +  2ig(aJfOpf  +  aJbOpb)flif -Figflrfflpf  +  /&,  (B26) 

^  =  (-iCDoi  -  rti/2)aib  +  ip*aif  +  2/g(aJfapf  +  ajbapb)aib  +  iga^a^  +  i'Cib,  (B27) 

where  £jv  =  yT^jve“'c'U  +  yT^Wjv  (j  =  p,s,i  and  v  =  f,b). 

Similar  to  the  previous  case,  Eqs.  (B22)-(B27)  can  be  used  to  find  the  solution  of  signal  and 
idler  fields  inside  the  cavity  ajv  (J  =  5,  i  and  v  =  /,  b),  through  which  we  obtain  the  signal  and  idler 
fluxes  as 


(gam)) 


1  CS 1  ll 

2  71  J-o 


/-|-oo  r 

\Ml3(o)\2  +  \Ml4(o)\: 

-00 


rPir, 


eiA  ts 


r*+°° 


2  n  J-o 


|M3i(o))|2  +  |M32(©)p 


do, 

do, 


(B28) 

(B29) 
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and  the  pair  correlation  as 


Pc(ts,ti)  =  (fj (ts)fs(ts)fi(ti)}  -  {/}&)/*&)){/! (td/ifa)) 


=  resn 


i 

-  /  {rts  [Mn{(o)M^x((o)  +  Mi2(a>)M22((o)\  -  M^co)}  dco 

Z7t  J — oo 


where  is  the  elements  of  the  matrix  M(co)  which  is  given  by 


M~\(0 )  =  - 


Ua  -  rts/2 

0  ^ 

m 

ico  rts/2 

0 

-fcflpb 

ig>l  f)2 

0  i(0 

-  rti/2 

-W 

0 

ig'K  b)2 

-ipi 

ico  -  rti/2 / 

(B30) 


(B31) 


Equations  (B28)  and  (B29)  provide  the  theoretical  curves  of  individual  signal  and  idler  fluxes 
given  in  Fig.  4a-d  in  the  main  text.  Equation  (B30)  shows  that  the  photon-pair  correlation  is  only 
a  function  of  time  difference  T  =  ts  —  ti,  pc(hdi)  =  Pc( t).  Integrated  pc{ t)  over  r,  we  obtain  the 
photon-pair  emission  flux  Rc  =  pc(t)dt,  which  is  given  as  the  theoretical  curves  in  Fig.  4e-h 
in  the  main  text.  The  complete  expressions  of  Eq.  (B3)-(B5)  are  used  in  the  theoretical  calcula¬ 
tion,  where  the  effective  mode  volumes  and  the  mode  overlap  factors  are  obtained  from  the  finite 
element  simulation  of  the  device.  A  Kerr  nonlinear  coefficient,  02=  5.0  x  10-5  cm2 /GW,  is  used 
in  the  theoretical  calculation.  The  cavity  parameters,  such  as  Toj,  Tej,  and  /3 j  ( j  =  p,s,i ),  are 
measured  from  the  laser-scanned  transmission  spectrum  of  the  passive  cavity.  The  Pi  mode  at 
1584.1  nm  is  outside  the  scanning  range  of  our  lasers  and  its  cavity  parameters  cannot  be  obtained 
through  independent  measurement.  We  approximate  its  intrinsic  Q  and  mode  splitting  with  those 
of  the  nearest  Pi  mode  at  1564.9  nm  since  the  properties  of  cavity  modes  remain  similar  within 
each  mode  family.  Moreover,  we  approximate  its  external  coupling  with  that  of  the  P2  mode  at 
1553.6  nm  because  of  their  similar  external  couplings  as  observed  from  the  PL  spectra  shown  in 
Fig.  3b  and  c  of  the  main  text.  Figure  4d,h,i  of  the  main  text  show  that  these  parameters  provide  a 
good  approximation  for  the  cavity  properties  of  the  Pi  mode  at  1584.1  nm. 


Appendix  C:  Statistics  of  coincidence  photon  counting 

In  this  section,  we  provide  the  theoretical  description  of  the  statistics  of  coincidence  photon 
counting  with  gated  single  photon  detectors  (SPDs).  Assume  a  signal  and  idler  photon  strike  the 
detectors  at  time  ts  and  /,■  with  a  probability  density  of  pph(ts,ti)-  The  SPDs  then  produce  a  pair 
of  photoelectric  events  which  are  recorded  by  the  coincidence  counter  at  a  later  time  t\  =ts-\-  T\ 
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and  t2  =  tj  +  T2  (Fig.  8),  where  X\  and  T2  include  the  electron  transit  time  inside  the  SPDs,  the 
propagation  time  from  the  SPDs  to  the  coincidence  counter,  and  time  delays  introduced  by  the 
coincidence  counter,  with  a  joint  probability  density  of  pj(x\,  T2).  The  joint  probability  density  of 
the  photoelectric  events  is  thus  given  by 


PSlit\,t2 )  =  ?1§2  //  pph(txJl)pj(Xl.X2)dT\dT2 

/+°° 

Pphih-h  ~  X[) ) P\Rf(X[) )d T/), 

-OO 


(Cl) 


where  £,j  ( j  =1,2)  are  the  detection  efficiencies  of  the  two  SPDs,  and  we  have  used  the 
fact  that  Pph{tsJi)  =  pph(ti  —  ts)  which  is  only  a  function  of  the  time  difference  between  the 
signal  and  idler  arrival.  For  a  correlated  photon  pair,  Eq.  (B30)  indicates  that  pph(tsdi )  = 
as OLj ( fj ( ti )  ( ts ) fs ( ts )  fi (ti))  where  as  and  a,  are  the  propagation  losses  of  signal  and  idler  pho¬ 

tons  from  the  coupling  waveguide  (a  tapered  fiber  in  our  case)  to  the  SPDs.  In  Eq.  (Cl),  p\rf(x) 
is  the  instrument  response  function  (IRF)  of  the  coincidence  counting  system,  given  as 

/*+°° 

PIRf(^d)  =  /  Pj('tl,Xi+TD)dTl-  (C2) 

Jo 


The  IRF  is  independently  measured  by  launching  identical  ultrashort  pulses  onto  the  SPDs  and 
recording  the  coincidence  spectrum  (see  Method  in  the  main  text).  In  out  coincidence  counting 
system,  the  IRF  is  dominantly  induced  by  the  timing  jitters  of  the  two  SPDs.  As  the  detector  timing 
jitter  varies  with  its  quantum  efficiency  (QE),  the  IRF  is  recorded  at  different  QEs  for  specific 
experiments  (QE  is  15%  and  25%  for  Fig.  4  and  5  in  the  main  text,  respectively).  In  general,  the 
experimentally  recorded  IRF  is  well  described  by  a  Gaussian  function  since  the  timing  jitters  of 
the  SPDs  follow  a  Gaussian  distribution. 

As  a  detection  event  can  be  either  a  photoelectric  event  fired  by  a  photon  or  a  detector  dark 
count,  the  joint  probability  density,  p(t\ ,  ti),  that  the  two  detectors  produce  a  pair  of  events  at  time 
t\  and  tj  (Fig.  8)  is  thus  given  by 


p{h,h)  =  Psi{hdi)  +  PsD{hdi)  +  PDi{t\di)  +  PDDijxdi),  (C3) 


where  psd ,  Pdi,  and  poo  are  the  joint  probability  densities  of  an  event  pair  of  signal— dark  count, 
dark  count— idler,  and  dark  count— dark  count,  respectively,  produced  in  the  two  separate  SPDs. 
The  generation  of  a  detector  dark  count  is  independent  of  an  event  generated  in  another  detector 
(either  a  photoelectric  event  produced  by  an  incoming  photon  or  a  detector  dark  count).  As  a 
result,  psD(t\di)  =  Ps(h)pm(t2),  PDi{t\d2)  =  PDi(h)pi(t2),  and  pDDihdi)  =  PDi(h)PD2(t2), 
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FIG.  8:  a.  Schematic  of  coincidence  photon  counting.  A  signal  and  idler  photon  arrive  at  the  detectors  at  time  ts  and 
respectively,  with  a  probability  density  of  pph{ts,ti).  The  two  detectors  produce  a  pair  of  events  which  are  recorded 
by  the  coincidence  counter  at  time  t\  and  t2,  with  a  probability  density  of  p(t  1^2).  b.  The  time  relationship  of  the 
events  arriving  at  the  coincidence  counter.  Each  event  arrives  within  a  gating  window  centered  at  tjo  (7  =  1 , 2)  and 
with  a  width  of  Tq.  The  two  input  channels  of  the  coincidence  counter  has  a  time  delay  of  7^. 


where  ps(t),  pi(t ),  and  pDj(t )  (j  =1,2)  are  the  probability  densities  of  detection  events  from 
signal,  idler,  and  dark  counts,  respectively,  produced  in  the  two  detectors.  In  experiment,  ps(t)9 
pi(t ),  and  pDj(t )  are  equal  to  the  average  detection  rates  of  the  related  events,  which  can  be 
measured  independently.  Consequently,  Eq.  (C3)  becomes 

p(h,t2)=  Psi{t2-ti)+Ps(h)PD2(t2)+PDi{h)pi{t2)+PDi{h)pm{t2)-  (C4) 

Therefore,  the  coincidence  at  x  =  t2  — 1\  within  a  small  delay-time  interval  of  dx  is  given  by 

rho+Ta/2  rtio+TG/2 

p(x)dx  =  /  dt\  /  dt2p(t\.t2) 

Jt\Q-TG/2  Jt2o-Ta/2 

=  di{TG-  |t-  Td\)  [psi(?  ~  Td)  +  PsPD2  +  PdiPi  +  PD1PD2} ,  (C5) 

for  \z-Td\  <  Tq  and  p  ( r )  =0  otherwise.  In  experiment,  dx  corresponds  to  the  time  resolution 
of  the  coincidence  counter,  which  is  4  ps  in  our  system.  t\o  and  t2Q  are  the  centers  of  the  gating 
windows  in  the  two  channels  (Fig.  8).  As  a  result,  the  delay-time-dependent  coincidence  counting 
histogram  (we  term  as  the  coincidence  spectrum  here  for  simplicity)  is  obtained  accordingly  as 

C(x)dx  =  RQTp(x)dx,  (C6) 

where  Rq  is  the  clock  frequency  of  the  SPDs  and  T  is  the  data  collection  time. 

Equations  (Cl),  (C5),  and  (C6)  provide  a  complete  description  of  coincidence  photon  counting, 
from  which  we  can  obtain  the  true  coincidence  spectrum  Ct(x )  and  the  accidence  coincidence 
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spectrum  Ca( t)  as 


Cj{z)dz  =  RqT dz  (TG  -  \z-Td\)  [psi{z  —  Td)  - psPi] 

/+o° 

Pc(z  —  Td  —  T£))pirf(t£>)Jt£),  (C7) 

-OO 

CA{z)dz  =  RoTdziTG  —  \z  —  Td\)  (ps  +  Pd\ )(pi  +  Pdi),  (C8) 


where  pc(z)  is  given  in  Eq.  (B30).  The  theoretical  coincidence-to-accidental  ratio  (CAR)  shown 
as  the  red  curves  in  Fig.  4e-h  of  the  main  text  is  obtained  as 


CAR  = 


(C9) 


where  Tfwhm  is  the  full  width  at  half  maximum  of  the  coincidence  spectrum.  In  Eqs.  (C7) 
and  (C8),  ps  =  ^as(fj(t)fs(t))  +  pRs  and  pi  =  ^20Ci(fJ (t)fi(t))  +  pRi,  where  and 

(fj ( t)fi(t )}  are  given  by  Eqs.  (B28)  and  (B29),  and  prs  and  pr-  are  the  probability  densities  of 
detection  events  from  Raman  noise  photons  at  the  signal  and  idler  channels,  respectively  (See 
the  next  section  for  the  details  of  characterizing  Raman  noises).  On  the  other  hand,  Eq.  (C8) 
shows  clearly  that  the  accidental  coincidence  can  be  accurately  calculated  from  independently 
measured  ps,  pi,  and  poj  (j  =  1,2).  Therefore,  using  these  measured  parameters  together  with  in¬ 
dependently  recorded  true  coincidence  we  can  reproduce  the  experimentally  recorded  CAR.  This 
theoretical  approach  explains  well  the  saturation  of  CAR  at  low  pair  fluxes  and  is  commonly  used 
in  literature31.  It  is  shown  as  the  black  dashed  curves  in  Fig.  4e-h  of  the  main  text. 

The  theoretical  curves  shown  in  Fig.  5  of  the  main  text  are  described  by  Eq.  (C6)  (normalized 
by  its  peak  value),  in  which  pph(ts,ti )  =  its) f s(ts)  fi{U))  obtained  from  Eq.  (B30)  and 

Pirf(z)  is  given  by  a  Gaussian  fitting  to  the  experimentally  recorded  IRF.  The  cavity  parameters 
used  in  (fj ( ti)fj (■ ts)fs(ts)fi(ti )},  such  as  Top  rej,  and  /3/  ( j  =  p,s,  i),  were  independently  measured 
from  the  laser-scanned  transmission  spectrum  of  the  passive  cavity. 


Appendix  D:  Characterization  of  Raman  noise  produced  in  delivery  fibers 

To  characterize  the  noise  for  each  photon  mode,  we  recorded  the  photon  flux  when  the  center 
wavelength  of  the  tunable  bandpass  filter  was  tuned  away  from  the  photon  mode  so  that  the  latter 
was  completely  blocked  (see  Fig.  9).  The  pumping  condition  remained  the  same  as  what  we  used 
to  characterize  the  photon  pairs,  with  the  tapered  fiber  touched  down  to  the  nanoforks  and  the 
pump  wavelength  on  the  cavity  resonance.  Such  an  experimental  condition  ensures  that  the  pump 
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FIG.  9:  Schematic  of  the  experimental  setup.  Details  can  be  found  in  Fig.  2  of  the  main  text. 


power  distribution  inside  the  delivery  fibers  remains  unchanged.  The  photon  flux  recorded  under 
this  condition  is  not  produced  from  the  device  since  the  photon  wavelength  does  not  coincide  with 
any  cavity  resonance. 

Figure  10  shows  an  example  of  recorded  photon  fluxes  for  Pl2,  where  the  solid  and  open  circles 
show  the  cases  when  the  filter  centers  were  tuned  onto  and  away  from  the  photon  mode  pair, 
respectively.  The  noise  photon  fluxes  exhibit  a  linear  dependence  on  the  pump  power.  Figure  10 
shows  clearly  that  the  observed  noise  is  the  dominant  noise  source  accompanied  with  the  photon 
pairs.  To  further  verify  that  the  observed  noise  does  not  come  from  the  device,  we  directly  lifted 
the  taper  up  so  that  the  tapered  fiber  did  not  couple  to  the  device  and  the  pump  wave  propagated 
only  inside  the  delivery  fibers  (as  well  as  the  optical  components  in  the  experimental  system.  See 
Fig.  9).  Kept  other  conditions  unchanged  and  recorded  again  the  photon  fluxes,  we  obtained  the 
open  triangles  shown  in  Fig.  10,  which  are  very  similar  to  the  open  circles.  This  clearly  verifies 
that  the  noise  photons  are  not  generated  from  the  device  but  produced  inside  the  experimental 
testing  system. 

In  our  testing  system  (Fig.  9),  the  CWDM  MUX  has  a  band  isolation  greater  than  120  dB. 
As  a  result,  any  broadband  noise  produced  before  the  MUX  ( e.g .,  in  fiber  0 )  will  be  completely 
blocked.  On  the  other  hand,  the  CWDM  DEMUX  separates  the  pump  wave  from  the  photon 
pairs  and  no  noise  will  be  produced  in  fiber  3  after  the  DEMUX.  Moreover,  when  the  pump  wave 
is  nearly  critically  coupled  into  the  device,  the  pump  power  in  fiber  2  after  the  device  is  small. 
Therefore,  the  noise  should  be  dominantly  produced  inside  fiber  1  before  the  device,  which  has  a 
rough  length  of  ~5  m.  The  open  triangles  in  Fig.  10  have  higher  magnitudes  than  the  open  circles 
simply  because,  when  the  taper  is  lifted  up,  the  pump  power  will  increase  in  fiber  2  and  will  start 
to  produce  noise. 

Figure  10  shows  that  the  idler  channel  at  the  Stokes  side  of  the  pump  has  a  larger  noise  photon 
flux  than  the  signal  channel  at  the  anti-Stokes  side.  This  feature  together  with  their  linear  power 
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FIG.  10:  Photon  fluxes  recorded  for  P?,,  with  detector  dark  counts  subtracted.  The  solid  and  open  circles  show  the 
photon  fluxes  when  the  center  wavelengths  of  the  bandpass  filters  are  tuned  onto  and  away  from  the  Pi,  mode  pair, 
respectively,  with  the  tapered  fiber  touched  down  to  the  nanoforks.  The  open  triangles  are  recorded  when  the  tapered 
fiber  is  lifted  up  and  thus  does  not  couple  to  the  device.  The  blue  and  red  show  the  signal  at  shorter  wavelength  and 
idler  at  longer  wavelength,  respectively.  The  dashed  lines  are  only  for  eye  guidance. 

dependence  implies  that  the  noise  is  likely  to  be  the  Raman  noise  produced  inside  fiber  1.  To 
verify  this,  we  incremented  the  length  of  fiber  1  (with  the  tapered  fiber  not  coupled  to  the  device) 
and  monitored  the  noise  spectrum.  Figure  11a  shows  an  example  of  the  noise  spectra  recorded 
with  three  different  fiber  lengths  of  ~5  m,  ~29  m,  and  ~1  km.  The  similarity  of  the  spectra 
among  these  three  cases  indicates  clearly  that  they  share  an  exactly  same  physical  origin.  The 
noise  spectra  show  the  Raman  noise  characteristics  of  silica54,  with  the  Boson  peak  clearly  visible 
at  a  frequency  shift  of  ~  1.5  THz  away  from  the  pump. 

To  reveal  more  characteristics  of  the  noise  spectra,  we  tuned  the  pump  wavelength  to  the  left 
end  of  the  CWDM  and  recorded  again  the  noise  spectra  (Fig.  1  lb).  In  this  case,  the  noise  spectra 
primarily  show  the  Stokes  side.  Both  the  Boson  peak  and  Raman  gain  peak  are  clearly  visible,  with 
a  frequency  shift  of  about  1.5  and  13  THz,  respectively,  away  from  the  pump  frequency.  Tuning 
the  pump  wavelength  to  1570nm  we  recorded  the  noise  spectra  now  mostly  on  the  anti-Stokes 
side  (Fig.  11c).  Again,  we  observed  both  peaks  but  with  smaller  amplitudes.  These  observations 
clearly  verify  that  the  observed  noise  indeed  comes  from  the  spontaneous  Raman  scattering  inside 
the  delivery  fiber. 
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As  the  Raman  noise  is  very  broadband,  in  experiment,  we  used  a  pair  of  tunable  bandpass 
filters  with  a  3-dB  bandwidth  of  1.2  nm  in  front  of  the  SPDs  (Fig.  9)  to  cut  the  Raman  noise.  The 
Raman  noise  fluxes  associated  with  each  mode  pair  were  measured  with  the  filter  centers  tuned 
away  from  the  mode  pair  far  enough  that  the  latter  were  completely  blocked.  The  recorded  Raman 
noise  fluxes  were  used  in  Eq.  (C8)  for  finding  the  theoretical  CAR,  as  discussed  in  the  previous 
section.  Due  to  the  strong  frequency  dependence  of  Raman  noise  spectrum,  a  certain  uncertainty 
may  be  accompanied  with  the  recorded  Raman  noise  fluxes.  Note  that  the  Raman  noise  can  be 
further  suppressed  in  the  future  by  using  narrower  filters  since  the  correlated  photon  pairs  have 
much  narrower  spectra.  The  noise  can  also  be  reduced  by  shortening  the  length  of  fiber  1.  Such 
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FIG.  11:  a.  Raman  noise  spectra  recorded  with  different  lengths  of  fiber  7,  with  an  input  pump  power  of  about 
270  /I W.  The  spectra  were  recorded  at  different  transmission  ports  of  the  CWDM  DEMUX  (indicated  by  color),  for 
easy  suppression  of  the  pump.  The  solid  curves  and  dotted  curves  show  the  cases  with  a  fiber  length  of  ~29  m  and 
km,  respectively.  The  amplitudes  of  the  dotted  curves  were  reduced  by  a  factor  of  10,  for  easy  comparison  with 
the  thin  solid  curves.  The  inset  shows  the  case  with  a  fiber  length  of  ^5  m,  same  as  what  we  used  to  characterize 
correlated  photon  pairs.  The  pump  power  and  wavelength  are  the  same  for  these  three  cases.  The  spectral  cut-off 
beyond  ^1590nm  is  due  to  the  quantum-efficiency  cut-off  of  the  detector  used  in  our  sepctrometer.  b.  Raman  noise 
spectra  with  the  pump  located  at  1470nm,  showing  the  Stokes  side,  fiber  1  has  a  length  of  ^1  km.  c.  Raman  noise 
spectra  with  the  pump  located  at  1570nm,  showing  mostly  the  anti-Stokes  side.  Fiber  1  has  a  length  of  ~1  km. 
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Raman  noise  will  be  completely  absent  if  an  on-chip  laser39  is  used  as  the  pump. 
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